AFRL-MN-EG-TR-200 1-7076 


SURFACE  LAYER  THERMODYNAMICS  OF  STEEL 
PENETRATORS  AT  HIGH  AND  VERY  HIGH  SLIDING 
VELOCITIES 


JANUSZ  R.  KLEPACZKO 

UNIVERSITY  OF  FLORIDA 

GRADUATE  ENGINEERING  AND  RESEARCH 

CENTER 

1350  N.  POQUITO  ROAD 
SHALIMAR,  FL  32579 


JUNE  2001 


FINAL  REPORT  FOR  PERIOD  30  APRIL  2001  -  20  JUNE  2001 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  unlimited. 


20030225  087 


AIR  FORCE  RESEARCH  LABORATORY,  MUNITIONS  DIRECTORATE 

Air  Force  Materiel  Command  ■  United  States  Air  Force  ■  Eglin  Air  Force  Base 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other  than  in 
connection  with  a  definitely  Government-related  procurement,  the  United  States  Government  incurs  no 
responsibility  or  any  obligation  whatsoever.  The  fact  that  the  Government  may  have  formulated  or  in  any 
way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication,  or 
otherwise  in  any  manner  construed,  as  licensing  the  holder,  or  any  other  person  or  corporation;  or  as 
conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way 
be  related  thereto. 

This  technical  report  is  releasable  to  the  National  Technical  Information  Services  (NTIS).  At  NTIS  it  will 
be  available  to  the  general  public,  including  foreign  nations. 

This  report  has  been  reviewed  for  technical  content  and  assessed  to  be  accurate.  Due  to  extenuating 
circumstances,  it  is  not  practical  to  edit  this  report  for  format  and  it  is,  therefore,  approved  “as  submitted” 
by  the  author. 


FOR  THE  COMMANDER 


PaJL  fc-  -L.L  L 


Mr.  Paul  K.  Laird,  Technical  Director 
Ordnance  Division 


Dr.  Mary  L.  Hughes/Program  Manager 


Damage  Mechanisrrfs  Branch 


MfClGrk  E.  Herzog,  t^^al  Advisor 


Anyone  having  need  of  a  copy  of  this  report  should  first  contact  the  Defense  Technical  Information  Center 
(DTIC)  at  the  address  shown  below.  If  you  are  a  registered  DTIC  User  and  qualify  as  a  recipient  of  this 
document,  DTIC  can  provide  you  with  a  copy.  If  you  are  a  registered  DTIC  User  and  do  not  qualify  as  a 
recipient,  DTIC  can  submit  a  request  for  release,  on  your  behalf,  to  the  controlling  DoD  agency  for  their 
review  and  determination.  Please  do  not  request  copies  from  the  Air  Force  Research  Laboratory, 
Munitions  Directorate.  Requests  for  additional  copies  should  be  directed  to: 


Defense  Technical  Information  Center  (DTIC) 

8725  John  J.  Kingman  Road,  Ste  0944 
Ft  Belvoir,  VA  22060-62 1 8 

This  report  is  published  in  the  interest  of  the  scientific  and  technical  information  exchange.  Publication  of 
this  report  does  not  constitute  approval  or  disapproval  of  the  ideas  or  findings.  Do  not  return  copies  of  this 
report  unless  contractual  obligations  or  notice  on  a  specific  document  requires  its  return.  If  you  no  longer 
have  a  need  to  retain  this  document,  please  refer  to  the  Destruction  Notice  on  the  cover  page  for 
instruction. 


If  your  address  has  changed,  and  if  you  wish  to  be  removed  from  our  mailing  list,  or  if  your  organization  no 
longer  employs  the  addressee,  please  notify  AFRL/MNMW,  101  W.  Eglin  Blvd.,  Ste.  135,  Eglin  AFB,  FL 
32542-6810,  to  help  up  maintain  a  current  mailing  list. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


1 .  AGENCY  USE  ONLY  (Leave blank)  2.  REPORT  DATE  IT  REPORT  TYPE  AND  DATES  COVERED - 

—TITLE  AND  SUBTITLE  - ^ - JUNE  2001 -  FINAL  30  APR  01  -  20  JUN  01 

n;™D  SUBTITLE  5.  FUNDING  NUMBERS 

SURFACE  LAYER  THERMOMECHANICS  OF  STEEL  PENETRATORS  AT  HIGH  Contract  #: 

AND  VERY  HIGH  SLIDING  VELOCITIES  J0N: 

_  DC. 


6.  AUTHOR(S)  - - - - - 

JANUSZ  R.  KLEPACZKO 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  “  - - 

UNIVERSITY  OF  FLORIDA 

GRADUATE  ENGINEERING  AND  RESEARCH  CENTER 
1350  N.  POQUITO  ROAD 
SHALIMAR,  FL  32579 
('850')  883-9350 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AIR  FORCE  RESEARCH  LABORATORY,  MUNITIONS  DIRECTORATE, 
DAMAGE  MECHANISMS  BRANCH  (AFRL/MNMW) 

101  W.  EGLIN  BLVD.,  STE.  135 
EGLIN  AFB,  FL  32542-6810 

PROGRAM  MANAGER:  DR,  MARY  L.  HUGHES  £8501  882-7998 

11.  SUPPLEMENTARY  NOTES  - —  - 2 - 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT  - - - - 

DISTRIBUTION  A:  Approved  for  public  release;  distribution  unlimited. 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORINGfMONITORING 
AGENCY  REPORT  NUMBER 

AFRL-MN-EG-TR-2000-7076 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words)  •  “  - - * - - - — - - 

One  of  the  most  difficult  problems  in  the  mechanics  of  kinetic  energy  projectiles  penetrating  geological  or  cementitious  target 
materials  is  a  proper  understanding  of  the  frictional  properties  between  the  projectile  surface  and  the  target.  The  model 
developed  in  this  study  is  limited  to  a  simple  definition  of  the  state  of  the  projectile  surface  as  a  set  of  uniformly  distributed 
micro-asperities  with  the  same  active  height,  h  =  20  microns,  which  are  subject  to  fast  adiabatic  shearing.  Although  the 
heat  conduction  equation  has  not  been  solved  numerically  in  this  study,  a  useful  approximation  of  the  evolution  of 
temperature  in  the  bulk  material  was  applied.  This  approximation  permitted  finding  the  closed  form  solution  for  the 
homologous  bulk  temperature,  at  different  sliding  conditions.  Because  evolution  of  the  bulk  temperature  could  be  estimated, 
it  also  permitted  finding  the  evolution  of  the  coefficient  of  friction  as  a  function  of  the  sliding  velocity.  For  every  definition 
of  the  coefficient  of  friction  a  substantial  decrease  of  the  resistance  to  sliding  at  increasing  velocities  has  been  found.  An 
open  question  remains  as  to  the  role  of  the  hydrostatic  pressure.  For  example,  the  experimental  data  suggest  that  at  constant 
velocity  the  coefficient  of  friction  will  diminish  as  a  function  of  pressure  as  predicted  by  Eq.  (51).  Thus,  an  advantage  of  the 
model  presented  here  is  that  hydrostatic  pressure  can  be  taken  into  consideration  in  a  correct  way.  Albeit  implementation  of 
the  friction  model  in  the  form  of  Eq.  (51)  into  numerical  codes  is  rather  too  early,  preliminary  trials  should  be 
recommended. 


14.  SUBJECT  TERMS  - - " - - - 

Thermomechanics;  Penetrators;  Friction;  Wear;  High  Velocity;  Penetration;  Geological 
Targets,  Concrete;  Projectiles;  Coefficient  of  Friction;  Thermodynamics;  Asperities 


17.  SECURITY  CLASSIFICATION 

OF  REPORT 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

UNCLASSIFIED 

15.  NUMBER  OF  PAGES 

73 

16.  PRICE  CODE 


19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT 

OF  ABSTRACT 

UNCLASSIFIED 

UL 

jjianuaiu  i  Jim  400  intv.  4.-00  ICUJ 

Prescribed  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro.  WHS/DIOR.  Oct  94 


SURFACE  LAYER  THERMOMECHANICS  OF  STEEL  PENETRATORS 
AT  HIGH  AND  VERY  HIGH  SLIDING  VELOCITIES 

PREPARED  BY 

JANUSZ  R.  KLEPACZKO* 


University  of  Florida 

Graduate  Engineering  &  Research  Center 

1350  N.  Poquito  Road 
SHALIMAR  FL  32579 

April  2001 


*Visiting  Professor 

Laboratory  of  Physics  and  Mechanics  of  Materials 
Metz  University 

lie  du  Saulcy,  F-57045  Metz,  France 


Abstract 


One  of  the  most  difficult  problems  in  the  mechanics  of  KE  projectiles  penetrating 
geological  or  cementitious  target  materials  is  a  proper  understanding  of  the  frictional 
properties  between  the  projectile  surface  and  the  target.  The  model  developed  in  this  study  is 
limited  to  a  simple  definition  of  the  state  of  the  projectile  surface  as  a  set  of  uniformly 
distributed  micro-asperities  with  the  same  active  height  h  =  20  pm,  which  are  subject  to  fast 
adiabatic  shearing.  Although  the  heat  conduction  equation  has  not  been  solved  numerically  in 
this  study,  a  Useful  approximation  of  the  evolution  of  temperature  in  the  bulk  material  was 
applied.  This  approximation  permitted  finding  the  closed  form  solution  for  0b>  the 
homologous  bulk  temperature,  at  different  sliding  conditions.  Because  evolution  of  the  bulk 
temperature  could  be  estimated,  it  also  permitted  finding  the  evolution  of  the  coefficient  of 
friction  as  a  function  of  the  sliding  velocity.  For  every  definition  of  the  coefficient  of  friction  a 
substantial  decrease  of  the  resistance  to  sliding  at  increasing  velocities  has  been  found.  An 
open  question  remains  as  to  the  role  of  the  hydrostatic  pressure.  For  example,  the 
experimental  data  suggests  that  at  constant  velocity  the  coefficient  of  friction  will  diminish  as 
a  function  of  pressure  as  predicted  by  Eq.(51).  Thus,  an  advantage  of  the  model  presented  here 
is  that  the  hydrostatic  pressure  can  be  taken  into  consideration  in  a  correct  way.  Albeit 
implementation  of  the  friction  model  in  the  form  of  Eq.(51)  into  numerical  codes  is  rather  too 
early,  preliminary  trials  should  be  recommended. 
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1.  Introduction 


The  most  difficult  problem  in  the  penetration  mechanics  of  steel  projectiles  into  a  concrete 
target  is  a  proper  understanding  of  frictional  properties  between  the  projectile  surface  and  the 
perforated  medium.  Since  velocities  of  penetration  are  very  high,  of  the  order  of  1,500  m/s, 
and  the  pressure  on  the  projectile  surface  is  very  high  as  well,  this  unique  combination  of 
parameters  makes  the  problem  difficult.  The  initial  kinetic  energy  of  a  projectile  might  be 
very  high,  since  it  increases  proportionally  with  the  second  power  of  velocity.  Thus,  there  is 
no  linear  relation  between  those  two  quantities.  The  main  factor  that  governs  the  penetration 
distance  is  the  force  due  to  the  target  resistance  superimposed  on  the  friction  conditions 
present  between  the  projectile  and  the  target.  Many  physical  processes  may  occur  during  a 
high-speed  penetration;  for  example,  erosion  of  the  projectile  surface  accompanied  by  high 
temperatures  and  possible  phase  transformation,  on  or  near  the  steel  surface,  transportation  of 
the  melted  material  onto  the  shank,  injection  of  the  melted  material  into  the  target  medium 
and  chemical  reactions  on  the  surface  like  fast  oxidation.  Those  physical  processes  which  are 
involved  in  penetration  mechanics  cause  the  estimation  of  a  “mean”  coefficient  of  friction  in 
the  Coulomb  sense,  that  is  Ft  =  p  Fn  (where  Ft,  F„  and  p  are  respectively  the  tangential  force, 
the  normal  force  and  the  coefficient  of  friction)  to  become  extremely  difficult.  It  can  be  shown 
that  the  coefficient  of  friction  defined  in  such  a  manner  is  not  a  constant  but  depends  on  many 
factors  like  sliding  velocity,  normal  pressure,  composition  of  penetrated  material  (aggregate 
size  and  paste),  and  on  the  yield  stress,  rate  of  thermal  softening,  thermal  conductivity, 
specific  heat  and  melting  temperature  of  the  projectile  material.  In  conclusion,  a  more  detailed 
study  of  friction  conditions  on  the  projectile  interface  is  of  great  importance  to  derive 
alternative  friction  laws  based  on  the  principles  of  physics.  Recent  studies  on  abrasion  and 
penetration  (mostly  numerical)  [1,2,3],  clearly  indicate  the  prudent  direction  of  future  efforts. 

During  the  last  few  decades  substantial  progress  has  been  witnessed  in  understanding 
temperature-coupled  processes  of  plastic  deformation.  One  such  process  is  Adiabatic  Shear 
Banding  (ASB).  Many  publications  have  appeared  concerning  this  topic,  including 
monographs  -  see,  for  example  Reference  [4],  along  with  the  information  recorded  in 
References  [5-10].  The  mechanism  of  ASB  is  relatively  simple  -  the  work  of  plastic 
deformation,  converted  largely  into  heat,  softens  the  material,  leading  to  instability  and 
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localization  of  plastic  deformation.  The  process  is  repeated  in  a  loop  and  the  final  result  is  a 
very  narrow  band  of  severe  plastic  deformation  ranging  in  width  from  5  to  100  micrometers. 
As  it  has  been  shown  in  Reference  [10]  many  physical  factors  influence  triggering  and 
evolution  of  the  ASB.  A  complicated  interplay  between  strain  hardening,  thermal  softening, 
strain  rate  sensitivity,  thermal  conductivity  and  sometime  local  inertia  occurs  in  an  ASB.  In 
some  materials,  particularly  in  steels  and  titanium  alloys,  a  phase-like  transformation  can  take 
place  within  a  narrow  zone  of  the  highest  temperature.  The  ASBs,  with  or  without  phase 
transformation,  often  act  as  sites  of  Mode  II  fracture  initiation. 

It  has  also  been  shown  that  formation  of  ASBs  depends  on  the  initial  and  boundary 
conditions  [10,1 1];  in  particular,  the  velocity  imposed  on  a  deforming  body  is  one  of  the  main 
factors  influencing  the  evolution  of  an  ASB.  If  the  imposed  velocity  is  high,  on  the  order  of 
100  m/s,  the  Critical  Impact  Velocity  in  shear  is  reached  [12,13],  and  a  localization  of  plastic 
deformation  will  occur  on  the  surface  where  that  velocity  is  imposed.  Importance  of  the  roles 
ASBs  play  in  diverse  processes  like  rolling,  machining,  drawing,  impact  on  structures, 
ballistic  impact  and  fragmentation  of  metals  and  alloys  is  obvious.  One  recent  illustration  is  an 
analysis  of  the  coefficient  of  friction  in  high-speed  machining  [14].  Additionally,  it  has  been 
shown  that  application  of  the  ASB  concept  to  analysis  of  friction  at  velocities  from  1  m/s  to 
100  m/s  leads  to  reasonable  agreement  with  experimental  data  [15,16],  In  these  studies,  the 
adiabatic  shear  banding  approach  to  the  process  of  surface  sliding  is  extended  to  high  and  very 
high  velocities,  similar  to  those  encountered  in  penetration  mechanics. 

2.  Fundamental  assumptions  and  the  surface  algebra 

Since  every  surface  has  a  certain  level  of  roughness,  during  sliding  the  tops  of  the  highest 
asperities  will  be  loaded  by  a  high  shear  stress.  If  the  sliding  velocity  is  sufficiently  high  the 
tops  of  the  highest  asperities  will  be  cut  by  adiabatic  shear  banding.  The  debris  of  the  tops  is 
the  source  of  the  surface  erosion.  Since  the  localization  process  of  plastic  deformation  in  ASB 
produces  high  temperature,  the  active  asperities  will  become  the  “hot  points”  and  the  heat 
generated  will  be  transmitted  into  the  bulk  of  the  material.  Because  the  process  of  adiabatic 
cutting  will  be  repeated  with  a  high  frequency  at  random  places  on  the  surface,  the  bulk 
material  will  be  constantly  heated  until  the  sliding  is  ended.  This  mechanism  produces  very 


7 


high  temperature  gradients  close  the  surface  of  the  bulk  material,  developing  a  Heat  Affected 
Zone  (HAZ)  not  necessarily  associated  with  a  phase  transformation.  It  may  be  mentioned  that 
the  phase  transformations  in  metals  are  normally  associated  with  a  sufficiently  high  cooling 
rate.  The  origin  of  the  adiabatic  cutting  is  not  specified  here.  It  is  understood  that  in  the  case 
of  geologic  materials  like  concrete  the  mineral  particles  will  erode  a  projectile  surface  (steel, 
for  example). 

The  main  mechanism  assumed  here  which  leads  to  evaluation  of  the  surface 
thermodynamics  and  erosion  is  based  on  dynamic  plasticity  with  thermal  coupling.  Although 
the  plasticity-based  models  of  friction  are  not  new  [15,16,17],  the  concept  of  applying 
adiabatic  shear  banding  as  a  principal  mechanism  in  friction  analysis  is  relatively  new  [14]. 

In  general,  the  surface  topography  can  be  represented  as  a  set  of  vertical  micro-asperities 
with  randomly  distributed  highs  hn  and  with  a  constant  square  cross  section.  Such  an 
idealization  is  shown  in  Fig.  1 .  In  this  approximation  three  levels  of  micro-asperity  heights  can 
be  recognized:  first  are  the  micro-asperities  of  the  highest  level  (red),  second  are  the  (green) 
micro-asperities  of  the  mean  height,  and  finally  there  is  the  subset  of  the  lowest  micro¬ 
asperities  (blue).  The  total  number  of  asperities  is  N  over  the  surface  S,  so  the  mean  density  of 
the  regular  array  of  micro-asperities  with  the  mean  square  cross  section  A2  is  n  =  N/S.  Since  S 
=  N  A  ,  the  density  of  all  micro-asperities  is  related  to  the  mean  side  dimension  A  by  the 
relation  n  =  1/A2. 

It  is  assumed  that  the  idealized  three  subsets  of  micro-asperities  have  different  functions 
during  the  sliding  process.  The  subset  with  the  highest  micro-asperities  is  cut  by  adiabatic 
shear  banding  and  at  the  same  time  it  carries,  to  some  extent,  the  normal  stress.  The  medium 
subset  supports  the  normal  stress  but  is  not  cut.  Finally,  the  subset  with  the  shortest  micro¬ 
asperities  rests  without  contact.  The  highest  micro-asperities  produce  the  “hot  spots”  shown  in 
red  on  the  top  view  of  the  surface  S  in  Fig.2. 

The  surface  topography  approximated  by  the  random  distribution  of  heights  must  involve  a 
statistical  analysis  of  the  adiabatic  cutting.  It  is  too  early  to  apply  such  an  approach,  and  in  the 
following  analysis  a  simplified  version  of  the  asperity  distribution  has  been  assumed,  that  is,  a 
uniformly  distributed  array,  Na,  of  the  active  micro-asperities  of  the  equal  height  h  which  fail 
by  adiabatic  shearing.  The  elementary  cell  (active  asperity)  with  the  active  micro-asperity  is 
shown  in  Fig.3.  In  this  approach  the  surface  S  is  covered  by  a  regular  array  of  square 
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elementary  cells  (active  asperities)  with  the  mean  dimension  A,  thus  S  =  Na  A2,  and  the 
density  of  such  defined  asperities  is  na  =  1/A2.  Every  elementary  cell  carries  one  vertical 
micro-asperity  of  the  cross  section  a2,  which  is  cut  by  ASB.  The  elementary  cell  surface  is 
only  partly  occupied  by  the  active  micro-asperity;  the  rest  of  the  surface  is  filled  by  the  lower- 
height  micro-asperities  invisible  in  this  approximation.  The  fractional  area  of  the  active  micro¬ 
asperities  can  be  defined  as  follows  [14], 


fa- 


fa}1 

\A  j 


fa*  1 


(1) 


In  general,  when  three  levels  of  height  are  assumed,  the  surface  S  is  divided  into  three 
domains: 

S  =  S0+Sa+Sb  '  (2) 


where  the  So  domain  is  the  surface  with  no  contact,  Sa  and  Sb  are  respectively  the  surfaces 
occupied  by  the  active  (adiabatic  cutting)  and  passive  (Coulomb  friction)  micro-asperities. 
Assuming  a  constant  number,  N,  of  all  asperities,  the  relation  between  particular  fractions  is 
given  by 


^  f nO  f sO  f na  f sa  "b  f nb  f s 


sb 


(3) 


with 


f  = *Ll 

Jn0  N 


fsO 


/  \  2 
'  C  ' 


f  =—  /  = 

5  J  na  jy  ’  ^ sa 


/  \  2 
'  a 

UJ 


f  f  -fA) 

?  J  nb  A7-  »  J  sb  A 

N  \ A) 


where  b  and  c  are  the  side  dimensions  of  the  micro-asperities  with  Coulomb  friction  and  with 
no  contact,  respectively.  Equation  (3)  can  be  rewritten  in  the  form 


1  = 


\AaJ 


(  u  ^2  f 


+ 


\AbJ 


+ 


\AcJ 


(4) 


where  Aa,AbandAc  are,  respectively,  the  mean  distances  between  particular  micro¬ 
asperities,  or 

1  =  naa2  +nbb2  +ncc2  (5) 

where  na,nb,nc  are  the  densities  of  the  respective  micro-asperities..  The  surface  fractions 
f sa  ^  fsb  are  *e  geometry  parameters  of  the  problem.  Assuming  further  that  a  =  b  =  c,  that 
is>  fsa  =  fsb  =  /,o  ’  Ecf(3)  becomes 
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1  =  /M(/«0  +fna+fnb)  (6) 

This  last  relation  is  important  since  the  number  fraction  of  the  particular  contacts,  that  is  the 
active  and  Coulomb  micro-asperities,  changes  as  a  function  of  pressure.  An  increase  of 
pressure  will  increase  the  number  of  active  micro-asperities  while  simultaneously  diminishing 
those  of  the  Coulomb-type.  Introducing  the  principle  of  mixture  and  the  influence  function,  p, 
of  pressure  (denoted  as  p(p)),  Eq.(6)  becomes 

i  =  fJp(P)f„+^-p(pWM  +/J]  (7) 

It  is  obvious  that  p(p)  must  be  an  increasing  function  of  pressure.  It  is  assumed  in  the 

simplified  modeling  that  the  effect  of  pressure  is  introduced  through  modification  of  Eq.(l)  as 
follows: 


Aa 


^  a  ' 


\AaJ 


fnaP(P) 


(8) 


where  fna  p(p)  _  ~s~-  ,  na/n  is  the  density  fraction  of  the  active  (adiabatic)  micro- 

aspenties,  0  <  fna  p(p)  <  1 ,  n  is  the  total  density  of  asperities.  In  this  approach  the  area 

fraction  remains  constant  but  the  density  of  the  active  micro-asperities  should  increase  as  a 
function  of  pressure.  At  the  same  time  the  mean  distance  Aa  between  the  active  micro¬ 
asperities  decreases.  The  evolution  of  the  density  na  is  assumed  in  the  form  of  the  following 
differential  equation 

dn. 


drM~Bn ■ 


(9) 


This  is  a  standard  evolution  equation  when  the  neighbor  entities  slow  down  the  growth  rate  of 
population;  M  is  the  constant  characterizing  the  ability  to  growth  of  the  active  micro-asperities 
and  B  is  the  constant  related  to  the  probability  of  interaction  between  neighboring  micro- 
asperities.  The  general  solution  of  Eq.(9)  is  given  by 


p  =  ~\n(M-Bna)  +  C  (10) 

The  constant  C  can  be  found  from  the  initial  conditions  na  =  na0  when  p-  p0.  The  final 
form  of  the  solution  is 
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_  M 

f  M 

~  B 

\~B~ 

Limit  of  the  solution  is 


Finally 


lim  na  ,  M 

- =  n  and  n  =  — 

p  ->oo  B 

=  1  ■ ~  (1  ■ - — )  exp(-  B(j>  -  p0 )) 

1 1  Yi 


The  solution  (1 3)  can  be  introduced  into  Eq.(8),  the  final  result  is 

Ml-  —  )exp(-5(/>-/?0))  (14) 

1Ya  \1Va  J  n 

Equation  (14)  will  be  used  later,  in  the  different  stages  of  this  study,  to  calculate  variety  of 
parameters. 


3.  Coefficient  of  friction 


By  definition  the  coefficient  of  friction  is  the  ratio  of  the  mean  tangential  force  Fx  to  the 
mean  normal  force  Fy  ,  thus 

F. 

Fy 

(15) 

Since  two  types  of  the  micro-asperities  in  contact  are  assumed, 
rewritten  as 

_  Fxa  +  Fxb 

“a”  and  “b”,  Eq.(15)  can  be 

Fya+F*. 

Introducing  tangential  and  normal  stresses  the  mean  forces  are 

(16) 

F,=N„a2ra+Nbb2Tb 

(17) 

Fy  -  Na  a2  tja  +  Nb  b2  crb 

(18) 

where  ta ,  tb  ,  ca ,  o  b  are  shear  and  normal  stresses  on  particular  micro-asperities.  Introducing 
Eqs.(17)  and  (18)  into  Eq,(15)  together  with  the  following  approximations:  a  ~  b, 
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^ b  Mb  ^  the  local  Coulomb  friction  and  p  —  <ja  ~  o~b .  that  is  the  normal  pressure 

the  generalized  coefficient  of  friction  is  given  by 

”a(P)*a  ,  nh  (P) 

09) 

As  expected,  the  total  coefficient  of  friction  in  this  model  combines  two  terms,  the  first  one  is 
due  to  shear  stress  understood  as  an  result  of  plastic  deformation  of  the  active  micro¬ 
asperities,  isothermal  or  adiabatic,  the  second  one  is  due  to  standard  Coulomb  friction 
occurring  in  non-deformed  plastically  asperities.  Because  plasticity  is  rate-dependent,  the 
shear  stress  and  the  generalized  coefficient  of  friction  must  be  a  function  of  sliding  velocity  V. 
It  is  interesting  to  note  that  both  terms  depend  on  the  normal  pressure  p.  Since  the  population 
of  the  active  micro-asperities  increases  as  a  function  of  pressure,  as  it  is  predicted  by  Eq.(13), 
and  at  the  same  time  the  Coulomb  micro-asperities  diminishes  with  the  same  proportion  (total 
n  assumed  constant),  in  the  limit,  plasticity  will  dominate  at  high  pressures.  An  inverse 
situation  is  expected  for  low  normal  pressure,  the  Coulomb  friction  will  dominate.  The  next 
step  in  simplifications  is  to  neglect  the  Coulomb  friction  and  to  analyze  the  plasticity  term. 
Since  n  =  na+nb  Eq.(19)  can  be  written  as 


M 


_  na(p)Tg 
pn 


(20) 


Introduction  of  Eq.(13)  into  Eq.(20)  yields  the  coefficient  of  friction  based  on  plasticity 

T_ 


M  = 


_  "a 


p  ^-0“ exp(-£ (p - p0 ) 


(21) 


However,  in  this  approximation  the  shear  stress  is  assumed  as  sliding  resistance  without 
specification  of  the  fractional  area  of  the  elementary  cell.  The  effective  shear  stress  teff  should 
be  reduced  with  the  following  proportion,  [13], 


=  V 


VA  J 


(22) 


Finally,  the  coefficient  of  friction  due  to  plasticity  is  given  by 


M  = 


^  [Ml -/a0)  exp  (-#  (/>-/><>)] 


(23) 
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with  fa0  =  na0  /  n  and  the  subscript  eff  is  omitted.  The  most  important  part  of  this  analysis 
is  the  ratio  x/p.  It  is  assumed  that  the  normal  pressure  p  is  a  fraction  x  of  the  yield  stress  in 
compression  ay,  thus  p-x<jy  and  Eq.(23)  becomes 

M  =  Y  l-(l-/a0)exp  -Bx0a  —  -1  (24) 

X°y\AJ  [  {  _  Vo  ))_ 

where  p0  =  x0  cr y.  In  order  to  show  the  effect  of  the  number  of  active  contact  points  on  p  the 

term  in  the  square  brackets  was  analyzed  with  the  following  initial  conditions:  fao  =  0.1  for  x0 
=  0.01  and  fmax  =  0.99  for  x  =  1,  B  =  5.04655*  10'3 ,  the  result  is  shown  in  Fig.4.  As  expected, 
at  high  pressure  the  saturation  level  of  the  number  of  active  asperities  is  reached. 

The  next  important  step  is  discussion  of  the  ratio  x  /  ay .  It  is  obvious  that  a  yield  condition 

must  be  introduced.  Many  possibilities  can  be  envisaged.  For  example,  in  [11]  the  Prager- 
Ziegler  kinematic  hardening  was  assumed  where  shear  stress  is  independent  of  the  normal 

stress.  However,  if  the  Huber-Mises  yield  condition  is  assumed,  that  is  xy  =  cry/  -\[3 ,  and 
the  shear  stress  is  assumed  to  increase  proportionally  with  the  normal  stress,  that  is  x-xxy  , 
this  results  in  a  Coulomb-type  equation,  where  the  coefficient  of  friction  is  given  by 

1  (a  VT  «—  fx  Y| 

M  =  ~j=  —  l-(l-/fl0)exp  -73 Bx0t  —1  (25) 

The  limit  of  the  coefficient  of  friction  is  reached  for  x  =  1  and  then  p  =  0.577  (for  a/A  =  1). 
The  result  of  calculations  for  this  model  of  friction,  with  xy  =  520  MPa  and  for  different 
fractional  area,  as  a  function  of  non-dimensional  pressure  is  shown  in  Fig.  5.  Another 
possibility,  a  very  probable  one,  is  that  the  active  micro-asperities  are  always  plastified,  that  is 

x  =  xy  and  pressure  increases  in  proportion  to  the  yield  stress,  that  is  p  =  J 3  xxy  .  This 

approach  to  the  coefficient  of  friction,  which  should  be  rather  called  friction  parameter, 
because  the  definition  the  coefficient  of  friction  in  the  Coulomb  sense  is  not  obeyed,  can  be 
written  as  follows 

1  r  yl”  f  f  \V 

M=  7=  ~  l-(l-/a0)exp  -V3 Bx0x  —  -1  (26) 

xV3(A;[_  ^  y0  )) 
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The  result  of  the  numerical  calculation  for  the  same  conditions  as  in  the  previous  case  is 
shown  in  Fig.6.  This  time  the  friction  parameter  reaches  relatively  high  values  and 
substantially  decreases  as  a  function  of  pressure. 

In  conclusion,  the  ratio  t/p  is  a  decisive  variable  in  definition  of  “friction”,  or  rather 
resistance  to  sliding,  in  all  problems  involving  plasticity,  visco-plasticity  and  temperature- 
coupled  visco-plasticity.  Experimental  results  collected  in  [15]  and  [16]  clearly  indicate  that 
the  resistance  to  sliding  substantially  diminishes  when  the  sliding  velocity  increases.  Further 
experimentation  is  of  great  importance,  especially  at  very  high  sliding  velocities,  ie.,  above 
1 00  m/s,  to  shed  some  light  on  dynamic  resistance  to  friction. 

4.  Estimation  of  strain  rates  at  high-velocity  sliding 

It  is  well  known  that  at  high  strain  rates  all  materials  show  very  high  rate  sensitivity.  The 
range  of  strain  rates  expected  in  a  particular  problem  indicates  a  level  of  approximation 
needed  in  constitutive  modeling.  Since  the  heights  h  of  the  micro-asperities  under 
consideration  may  be  of  an  order  in  between  few  micrometers  to  say  few  tens  of  a  millimeter 
and  the  whole  range  velocities  of  interest,  say  from  1.0  m/s  to  about  1500  m/s,  the  range  of 
strain  rates  is  very  large.  In  order  to  estimate  this  range  a  simple  calculation  was  performed, 
the  results  are  shown  in  Fig.7  and  Fig.8.  The  shear  strain  rate  defined  by  standard  relation 

r  V 

r=J  ^ 

The  range  of  asperity  heights  is  assumed  as  5.0  pm  to  100  pm  (0.1  mm).  In  Fig.7  the  velocity 
was  assumed  as  low  as  1.0  m/s  and  the  range  of  strain  rates  is  in  between  104 1/s  (h  =  0.1  mm) 
and  106  1/s.  The  strain  rates  are  surprisingly  high  even  for  sliding  velocity  as  low  as  1.0  m/s. 
Fig.8  shows  the  logarithm  of  strain  rate  in  the  of  velocity  range  200  m/s  <  V  <  1600  m/s  as 
function  of  micro-asperity  heights.  Again,  strain  rates  are  extremely  high,  from  about  105  1/s 
(h  =  0.1  mm)  to  about  108  1/s  for  smaller  h. 

The  main  conclusion  reached  on  the  basis  of  those  calculations  can  be  formulated  in  the 
following  statement:  constitutive  relations  using  standard  dynamic  plasticity  cannot  be 
applied  to  find  the  resistance  to  sliding  friction  at  very  high  strain  rates.  Alternatively, 
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constitutive  relations  based  on  pseudo-viscosity  using  electron  or  phonon  viscosity  and 
supersonic  dislocations,  or  even  hydrodynamic  modeling,  might  prove  useful. 

5.  Constitutive  modeling  at  very  high  strain  rates  -  pseudo-viscosity 

It  is  commonly  assumed  in  constitutive  modeling  of  dynamic  plasticity  based  on  dislocation 
dynamics  that  at  specific  temperature  the  effective  shear  stress  required  to  activate  dislocation 
motion  without  assistance  of  thermal  vibration  of  atoms  is  the  limit  of  standard  visco¬ 
plasticity  approach.  This  level  of  shear  stress  is  called  the  mechanical  threshold,  [18].  This 
mechanical  threshold  is  reached  at  absolute  zero  temperature  or  at  very  high  specific  strain 

rate,  typically  ~5xl03  1/s,  then  at  T  =  0  K,  t  =  t0  and  f  =  fco .  Above  the  threshold  strain 

rate  rco  there  is  an  increase  in  defects  (dislocations)  together  with  a  substantial  increase  of 

the  mean  dislocation  velocity  [19].  High  dislocation  velocity  triggers  different  micro¬ 
mechanisms  of  plastic  flow  like  electron  or  phonon  viscosity,  which  leads  to  the  rate- 
dependent  threshold  stress,  [19].  In  other  words,  when  the  applied  stress  x  exceeds  the 
threshold  stress,  t  >  z ,  T  >  0,  the  dislocation  velocity  is  controlled  by  energy  dissipation  via 
electron  and  phonon  viscosities  when  dislocations  move  through  the  lattice.  It  is  commonly 
assumed  that  this  dissipation  is  linear  and  viscous  and  the  force  F  acting  on  a  mobile 
dislocation  of  length  lm  is  proportional  to  the  mean  dislocation  velocity  with  the 
proportionality  constant  B  called  the  viscous  drag,  such  relations  may  be  written  as 

F  =  Bv  F  =  zlm  b  lm=lljp^  (28) 

where  pm  and  B  are  respectively  the  mobile  dislocation  density  and  viscous  drag,  x  the  is  shear 
stress  on  a  dislocation,  b  is  the  magnitude  of  Burger’s  vector  (b  =  2.48xl0'8  cm  for  iron). 

The  Orowan  relation  relates  f ,  the  shear  strain  rate,  with  the  mobile  dislocation  density  and 
the  mean  dislocation  velocity 

r  =  bpmv  (29) 

The  shear  stress  can  be  written  as 
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B 


B 


(30) 


T0  + 


B  v 

Lb 


or  and 


T  =  T0  + 


b2  yJ~Pn 


-f  with  1)  = 


b2ip„ 


with  r),  the  pseudo-viscosity.  Finally,  at  strain  rates  higher  than  the  threshold  strain  rate,  and 
the  threshold  stress,  the  flow  stress  can  be  written  as 

t  =  t0(T)  +  tj(T)T  (31) 

The  threshold  stress  and  the  pseudo  viscosity  must  be  temperature  dependent,  where  T  is  the 
absolute  temperature. 

At  very  high  strain  rates  dislocations  approach  supersonic  velocity  C2,  C2  =  -^G/  ~p0  , 

where  G  is  the  shear  modulus  and  po  is  the  mass  density.  The  supersonic  approximation  to  the 
dislocation  drag  has  been  proposed  in  [1 9]  in  the  form 


F  = 


Bv 


-11/2 


and 


M^)2 


C. 


2  J 


r  =  r0  + 


BY 


b2  yfPn 


-1-1/2 


1- 


\b pm  C2  j 


(32) 


with  the  critical  strain  rate  f  cril  =  b  pm  C2 


V  =  Po 


r  r 

\b  Pm  ^2  J 


-1/2 


and  generalized  pseudo-viscosity  rj,  r  =  r0  +  r/Y 


with 


(33) 


It  is  interesting  to  analyze  theoretical  limits  of  the  viscous  drag  approximation.  Assuming  the 
elastic  wave  speed  in  shear  mode  €2  =  3.0  mm/ps  the  critical  shear  strain  rate  for  iron  may  be 
estimated  from  the  linear  relation  (  Eq.(29))  Ycril  =  7.44  *10'3  pm  ,  where  [pm]  =  [1/cm2].  For 

example,  for  the  mobile  dislocation  density  108  1/cm2  the  critical  strain  rate  is  7.44*  105  1/s 
and  for  1010  1/cm2  the  critical  strain  rate  reaches  value  of  7.44*  107  1/s.  These  estimations 
indicate  that  the  strain  rates  during  penetrator  sliding  are  very  close  to  the  limiting  strain  rates 
predicted  by  the  theory  of  dislocation  dynamics.  To  complete  this  part  of  the  analysis  Eq.(33) 
has  been  analyzed  for  different  values  of  the  mobile  dislocation  density  and  the  result  is 
shown  in  Fig.9. 

The  value  of  r|o  was  found  from  the  direct  impact  shear  experiments  on  mild  steel,  [19],  and 
earlier  experimental  results,  [21],  which  confirm  the  value  of  the  pseudo-viscosity  at  0  K,  rjo  = 
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2.65*10  Pa*s.  Theoretical  curves  show  that  the  viscous  approximation  with  the  constant 
value  of  r|o  is  correct  at  high  densities  of  mobile  dislocations,  say  above  108  1/cm2. 

Under  this  assumption  Eq.(3 1)  has  been  modified  to  account  for  the  temperature  effects 

Gf©) 

r(0,D  =  — -f0+7(®)r  t  >f0  (34) 

Go 

where  ©  =  T/Tm  is  the  homologous  temperature  and  Tm  is  the  melting  point.  Because  the 
thermal  softening  of  the  crystalline  lattice  diminishes  resistance  to  plastic  flow,  the 
mechanical  threshold  stress  at  0  =  0  (T  =  0)  is  normalized  by  the  relative  changes  of  the  shear 
modulus  of  elasticity  as  a  function  of  temperature,  Go  is  the  shear  modulus  at  0  =  0.  In  the 
literature  many  approximations  to  the  relative  changes  of  the  shear  modulus  as  a  function  of 
temperature  have  been  proposed.  Here,  the  simplest  approximation  that  fits  relatively  well 
experimental  data  is  introduced 

G(0)  , 

— —  =  1-0  0  < 0  <  1  (35) 
Go 

An  open  question  remains  with  regard  to  the  temperature  dependence  of  the  pseudo-viscosity. 
In  general,  the  pseudo-viscosity  must  be  a  decreasing  function  of  the  homologous  temperature 
0.  Since  the  surface  temperatures  of  KE  penetrators  rapidly  reach  the  melting  point,  an 
expression  for  the  pseudo- viscosity  must  also  take  into  account  the  viscosity  of  a  liquid  metal. 
The  following  phenomenological  expression  for  t](0)  satisfies  such  boundary  conditions  for 
temperature,  that  is  at  0  =  0,  t|(0)  =  r|0  and  at  0  =  1  ,  r|(l)  =  r|m 

Tl(  0)  1  r 

- =  — [l  -  tanh(yf  ©  -  5)]  (36) 

Vo  2 

where  A  and  B  are  constants.  Values  of  the  constants  that  satisfy  the  boundary  conditions  are 
A  =  4,  B  =  2.  The  relative  changes  of  the  pseudo-viscosity  as  a  function  of  the  homologous 
temperature  is  shown  in  Fig, 10.  The  model  is  based  on  the  ideal  plasticity  concept  (no  strain 
hardening)  and  both  relative  functions,  i.e.  the  shear  modulus  and  the  pseudo-viscosity,  are  a 
rapidly  decreasing  function  of  the  homologous  temperature. 


17 


6.  Thermal  coupling  in  adiabatic  conditions 


It  has  been  known  for  a  long  time  that  the  large  part  of  the  work  of  plastic  deformation  (on 
average  90%)  is  converted  into  heat.  In  the  19th  century  Tresca  remarked  on  “heat  lines” 
which  appear  during  forging.  An  adiabatic  process  of  plastic  deformation  occurs  when  the 
heat  conduction  between  neighboring  particles  and  surroundings  is  suppressed.  A  specific 
case  of  the  adiabatic  process  of  plastic  deformation  occurs  under  condition  of  normal  heat 
conduction,  but  at  high  strain  rates,  when  the  deformation  process  is  too  short  to  permit 
significant  heat  transfer  in  the  crystalline  lattice.  Adiabatic  deformation  triggers  a  specific 
instability  called  Adiabatic  Shear  Band,  discussed  in  the  section  1  of  this  study.  Many  papers 
have  been  published  on  occurrence  of  the  ASBs,  theoretical,  experimental  and  numerical.  It  is 
out  of  scope  of  this  paper  to  review  them.  However,  one  of  the  earliest  detailed  analysis  of 
adiabatic  process  of  deformation  and  analysis  of  the  stability  conditions  in  tension  can  be 
found  in  [22]. 

The  equation  of  heat  conduction  in  the  form  applicable  to  dynamic  plasticity  may  be  written 
as 


nr  5T  .  ar  ,d2T 
pc*  —  =  y#  HT,r  ,T)~ — x— — 

dt  dt  dy2 


where  y  is  the  direction  of  the  heat  conduction,  p,  Cv  and  X  are  respectively  the  mass  density, 
the  specific  heat  at  constant  volume  and  the  heat  conductivity  constant  (Fourier  constant),  P  is 
the  fraction  of  the  plastic  work  converted  into  heat,  p  «  0.9  .  In  case  of  X  -  0  the  process  is 
adiabatic.  When  the  homologous  temperature  is  introduced,  and  the  strain  hardening  is 
neglected,  the  energy  balance  equation,  Eq.(37),  is  transformed  into  the  following  form 


5©  p  X  d2® 

dt  pCJm  dt  pCv  dy2 


(38) 


Introducing  constitutive  relations  (34),  (35)  and  (36)  into  the  equation  of  heat  conduction  the 
estimation  of  the  Heat  Affected  Zone  (HAZ)  can  be  found  numerically.  In  case  of  the 
adiabatic  process  Eq.(38)  is  reduced  to  the  following  form,  [22], 


d®__  P 
PCvTm 


(39) 
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This  simple  ordinary  differential  equation  should  be  integrated  with  the  initial  conditions  0  = 
0o  for  T  =  0.  Since  the  simplified  constitutive  relations  do  not  depend  on  strain  the  separation 
of  variables  and  integration  is  direct.  Transition  from  the  isothermal  into  adiabatic  process  of 
deformation  was  analyzed  numerically  for  a  short  thin  tube  deformed  in  torsion  at  different 
rates.  For  steel  the  “50%”  transition  from  isothermal  to  adiabatic  deformation  was  found  to 
occur  at  f  «  40  Us .  No  doubt,  at  strain  rates  higher  than  100  1/s  the  process  of  deformation 
is  entirely  adiabatic. 


7.  Temperature  evolution  in  one  micro-asperity 


Assume,  as  the  first  approximation,  that  the  micro-asperity  is  in  adiabatic  state,  and  over  the 
height  h  shear  deformation  is  uniform.  Thus,  in  fact,  only  one  layer  of  equal  micro-asperities 
is  analyzed.  Assume  further  that  the  height  has  the  same  dimension  in  all  calculations,  h  =  20 
pm,  a  mean  value  of  ASB  thickness  observed  experimentally.  Widths  of  so  called 
“transformed”  ASB  varies  in  steel  from  ~5.0  pm  to  ~20.0  pm,  non-transformed  band 
thickness  varies  from  ~20.0  pm  to  even  100  pm.  The  scheme  assumed  in  numerical 
calculations  is  shown  in  Fig.  11  where  the  micro-asperity  is  deformed  in  shear  to  the  critical 
failure  strain  Tc.  The  point  of  shear  serves  as  the  source  of  a  fast  increase  of  temperature, 
usually  up  to  the  melting  level. 

It  is  of  interest  to  calculate  the  time  interval  it  takes  the  homologous  temperature  to  reach 
the  melting  point  or  failure  at  a  particular  velocity  of  sliding  V.  Equation  (39)  can  be  written 
as 


d®  = 


Vo 

2  t0  h 


(40) 


with 

/,(©)  =  l-tanh[2(2©-l)]  (41) 

where  f  =  V  /  h .  and  the  operator  (*}  introduced  into  Eq.(40)  means  that  if  ©2  >  1  then 


1-0  )  =  0  otherwise  1-02.  If  the  failure  criterion  in  shear  is  assumed  as  the  critical  shear 
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strain  rf,  the  time  interval  to  failure  is  given  by  tf=hTf/V.  Now  Eq.(40)  can  be  integrated 


and  the  interval  of  time  tj  found  for  a  specified  value  of  the  homologous  temperature  ©* 


©, 

h  =  J 

®0 


d@ 


Cl  t0 


1-©2W 


Vo 

2  t0  h 


4(0) 


v_ 
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P  4  Tm 


(42) 


Of  course,  the  most  important  parameter  of  the  problem  is  the  imposed  velocity  V.  The  set  of 
equations  has  been  analyzed  numerically  with  the  following  values  of  constants 
0o  =  0.1628  that  is  ©0  =  293  K/1800  K,  Tm  =  1800  K  for  iron 


|3  =  0.9 

p  =  7800  kg/m3 

Cv  =  500  J/kg  K  so  Ci  =  1.282052*1  O'10  m2/N 


h  =  20.0  pm 

x0  =  520  MPa  ,  this  value  of  the  mechanical  threshold  at  ©  =  0  is  estimated  from  the 
experimental  data  for  1018  steel  given  in  [20] 

T)0  =  2.65*  1 03  Pa*s,  this  value  of  the  pseudo-viscosity  at  0  =  0  is  estimated  from  the 
experimental  data  for  a  mild  steel,  [20,21],  r|RT  =  2.5*103  Pa*s 
The  results  of  calculations  are  shown  in  Figl  2,  Fig. 13  and  Fig. 14.  The  time  to  reach  particular 
homologous  temperature  ©  for  different  velocity  level  V  is  shown  in  Fig.  12a,  b  and  12c.Three 
velocity  ranges  were  calculated:  (a)  1.0  m/s<  Y  <10  m/s;  (b)  10  m/s  <  V<  100  m/s  and  (c) 
100  m/s  <  V  <  1000  m/s.  It  is  clear  that  the  melting  point  is  reached  very  quickly,  for  velocity 
10  m/s  the  melting  point  is  reached  after  about  40  ps,  for  V  =  100  m/s  after  about  0.8  ps  and 
for  V  =  1000  m/s  after  about  10  ns  ( 1  ns  =  10  9  s).  In  conclusion,  the  active  micro-asperities 
reach  almost  instantaneously  the  melting  temperature.  Different  view  of  the  numerical  results 
are  shown  in  Fig.  13  a,  b  and  c.  The  time  to  reach  particular  homologous  temperature  has  been 
calculated  for  three  different  ranges  of  velocities.  Again,  the  range  of  times  to  reach  particular 
homologous  temperature  diminishes  rapidly  with  velocity.  An  interesting  question  arises, 
what  is  reached  first  failure  or  melting?  The  mechanical  failure  criterion  defined  by 
tc  =  hYc  IV may  be  satisfied  before  the  melting  temperature  is  reached.  The  results  of  such 

analysis  are  shown  in  Fig  14  a,  b  and  c  in  the  form  of  time  to  reach  failure  criterion  with  the 
critical  shear  strain  Tc  =  2.0  or  the  melting  temperature.  At  low  and  medium-high  velocities 
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the  failure  dominates,  however  at  velocities  higher  than  about  200  m/s  (see  Fig.  14c)  the 
melting  point  is  reached  earlier  than  failure. 

8.  Heat  Affected  Zone  (HAZ) 


The  evolution  of  temperature  calculated  for  one  micro-asperity  is  one  of  the  heat  sources 
that  heats  the  bulk  material.  Other  heat  sources  may  be  in  the  form  of  micro-debris  rubbing 
the  surface  and  also  foreign  grains  of  perforated  material.  During  penetration  the  bulk  material 
is  constantly  heated  by  repeated  process  of  adiabatic  shearing  of  micro-asperities.  The  failure 
of  micro-asperities  is  repeated  with  a  high  frequency  and  every  new  failure  contributes  as  a 
“hot  point”  to  the  increase  of  temperature  of  the  bulk  material.  The  process  is  ended  after  time 
tp  =  Lp/V,  where  Lp  is  the  projectile  length.  It  is  interesting  to  calculate  the  frequency  of 
cutting.  Assuming  the  critical  strain  to  failure  Tc  =  2.0  and  the  micro-asperity  height  h  =  20 
pm  the  frequency  of  cutting  is  vc  =  V/  (h  Fc),  for  V  =  100  m/s  the  frequency  of  cutting  is  vQ  = 
2.5  MHz,  for  V  =  1000  m/s  the  frequency  reaches  vc  =  25  MHz.  With  this  frequency,  the  bulk 
surface  is  heated  practically  all  the  time  by  the  melt  on  the  surface.  In  order  to  estimate  the 
HAZ  the  heat  conduction  equation  should  be  solved,  Eq.(  38),  with  the  boundary  conditions 
©  =  1  on  the  bulk  surface  (y  =  0)  and  0  =  0O  inside  the  bulk  material  (for  steel  ©0  =  0.1628) 
the  initial  conditions  are  0  =  0O  for  t  =  0.  A  simplified  estimation  of  HAZ  has  been  proposed 
in  [14].  With  the  process  time  tp , which  is  the  characteristic  time  of  the  heat  conduction,  the 
HAZ  is  estimated  as  a  characteristic  thermal  diffusion  length  during  that  time.  The  result  is 


H  = 


a 


1  Lp 
pCv  V 


d/2 


(43) 


where  H  denotes  the  HAZ  and  a  is  the  characteristic  constants  of  the  problem  considered. 


Constant  a  has  been  discussed  in  [14]  and  estimated  as  a  =  1/J2x  on  the  basis  of  the 
solution  of  the  similar  boundary  value  problem  reported  in  [23],  It  is  interesting  to  note  that 
for  a  constant  projectile  length  the  thickness  of  the  HAZ  diminishes  in  proportion  to  1/  -Jv  . 
Assuming  thermal  conductivity  constant  X  =  54.0  W/m  K  (steel)  the  size  of  HAZ  has  been 
calculated  as  a  function  of  the  penetration  time  tp  =  L P/V  the  results  are  shown  in  Fig.  15.  The 
upper  part  shows  HAZ  in  millimeters  for  relatively  long  time  of  penetration,  the  lower  figure 
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for  the  short  times.  For  example,  when  Lp  is  assumed  500  mm  and  velocity  100  m/s  the 
penetration  time  is  relatively  long,  tp  =  5.0  milliseconds,  for  Lp  =  100  mm  and  V  =  1500  m/s, 
tp  =  66.7  microseconds.  The  estimation  of  HAZ  gives  the  dimension  range  from  about  200  pm 
to  about  5  mm.  This  estimation  is  not  directly  related  to  so  the  called  transformed  zone.  The 
transformed  zone  should  be  in  general  related  to  the  cooling  rate  of  a  material.  The 
transformed  zone  must  be  thinner  in  comparison  to  the  present  estimates  of  HAZ. 


9.  Estimation  of  the  bulk  material  temperature 


Since  the  HAZ  is  approximated  by  Eq.(43)  a  simple  method  can  be  used  to  approximate  the 
bulk  material  temperature,  [14].  The  energy  balance  between  the  input  energy  produced  by  the 
active  micro-asperities  and  the  absorbed  energy  by  the  HAZ  per  unit  time  and  over  the  surface 
S  leads  to  the  following  relation 


pCv{Tb-T,)SH  =  pT(V,T)V 


2 

^A  j 

l  n  j 
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(44) 


where  AT  =  Tb-T0  is  the  increment  of  temperature  due  to  the  heat  production  by  the  active 
micro-asperities.  Introducing  expression  for  H  (Eq.(43))  into  Eq.(44)  the  energy  balance 


becomes 
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(45) 


Equation  (45)  must  be  solved  for  0b  after  introduction  of  Eq.(14)  also  Eqs(34),(35)  and  (36). 
Analytic  solution  is  possible  only  for  a  constant  viscosity  p.  Denoting  the  following 
expressions  in  Eq.(45)  as 


D(fr.p.L.A  =  afltt(  L 
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the  simplified  notation  of  Eq.(45)  is  obtained 

©,=©0  +D(l-®2h+E) 

The  final  form  to  be  solved  is  in  the  quadratic  equation 


(46) 
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with  F  =  1  +  E . 

The  solution  of  Eq.(47)  is 


ri=const 


(47) 


(&  ^ 

Hi +  F 

[D  J 


=  0 


=^[±(1  +  4D(0O  +DF))'n  -l]  (48) 

The  solution  has  been  analyzed  within  the  lower  range  of  velocities  for  nine  fractional  areas 
and  fna  =  0.9.  The  results  are  shown  in  Fig.  16.  The  other  parameters  were  assumed  as  follows: 
ap  =  0.14324,  t0  =  520  MPa,  Lp  =  250  mm,  r|  =  2.5 *103  Pa*s,  h  =  20  pm.  The  other  constants 
are  given  in  the  previous  parts  of  this  paper.  The  explicit  solution  shows  a  very  rapid  increase 
of  temperature,  that  is  the  melting  point  is  reached  very  quickly.  The  main  reason  is  the 
constancy  of  the  pseudo-viscosity,  the  main  generator  of  heat.  Although  the  analytic  solution 
may  be  understood  as  the  upper  bound  of  temperature  increase,  a  complete  analysis  with  r|(@) 
can  provide  a  more  exact  temperature  evolution  in  the  bulk  material. 

The  results  of  such  a  complete  analysis,  which  involves  Eqs(14),(34),(35),(36)  and  Eq.(45) 
are  shown  in  Fig.  17.  The  roots  of  Eq.(45)  were  calculated  as  a  function  of  velocity  for  four 
ranges,  that  is  0  <  V  <  1.0  m/s,  Figl7a  ;  1.0  m/s  <  V  <  10  m/s,  Figl7b;  10  m/s  <  V  <  100  m/s, 
Fig.  17c  and  100  m/s  <  V  <  1000  m/s,  Figl7d.  A  very  consistent  result  has  been  found,  the 
homologous  temperatures  increase  steadily  as  a  function  of  velocity  up  to  different  levels  of 
saturation  depending  on  the  fractional  area  of  micro-asperities.  At  velocities  higher  than  10 
m/s  the  saturation  levels  of  the  homologous  temperatures  are  close  to  or  slightly  above  the 
melting  point  (0  =  1).  Within  the  velocity  range  100  m/s  <  V  <  1000  m/s  the  homologous 
temperatures  are  higher  than  the  melting  point  because  the  viscosity  itself,  this  time  real, 
produces  heat.  Of  course,  the  picture  is  of  qualitative  nature  since  a  true  viscosity  of  melted 
steel  is  unknown  at  present.  However,  the  trends  presented  in  Fig.  17  are  in  complete 
agreement  with  the  physical  intuition. 

Further  calculations  show  changes  of  the  bulk  temperature  as  a  function  of  the  fractional 
area  a/A  (asperity  ratio)  at  different  velocities.  The  results  are  shown  in  Figl8a,  b  and  c  for 
three  different  ranges  of  speeds.  The  fractional  area  has  an  important  influence  on  the 
temperature  evolution  in  the  bulk  material.  If  the  fractional  area  is  larger  the  heat  transfer  from 
the  micro-asperities  is  more  intense.  Again,  the  temperatures  of  the  melting  point  and  above 
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are  reached  very  quickly  when  the  fractional  area  is  increased.  At  very  high  velocities  the 
effect  of  the  fractional  area  is  reduced. 

10.  Calculation  of  friction 


The  fundamental  mathematical  formula  for  the  coefficient  of  friction  p  is  given  by  Eq.(23). 
By  substituting  for  r.(©,f )  from  Eqs(34),  (35)  and  (36)  an  explicit  formula  for  p  is  obtained 

=  <  '-(Q^))2  )  +  ~f-  [l-0-/,*,)exp(-i!Q>-p0)]  (49) 

where  the  operator  (*)  =  0  if®  >  1  otherwise  (*)  =  1  -  © \ .  The  solution  for  0b(V)  is  given 
by  Eq.(48)  and  Tj(®b)  by  Eq.(36). 

In  order  to  analyze  the  effect  of  sliding  velocity  on  the  coefficient  of  friction  Eq.(49)  was 
normalized  by  the  coefficient  of  friction  at  the  reference  velocity  V  =  1.0  m/s.  The  ratio  of  the 
coefficients  of  friction  is  given  by 


l-©^(F)  +  7(©)E/r0A 

l-®M(^)  +  <7o(®)r0/ro  * 


(50) 


where  Vo  1.0  m/s,  ®bo(Vo)  and  r|o(©)  denotes  respectively  the  evolution  of  the  homologous 
temperature  and  pseudo-viscosity  at  the  reference  velocity  V0  =  1.0  m/s.  The  ratio  has  been 


analyzed  numerically  and  the  results  are  shown  in  Fig.  19,  a,  b,  and  c  for  different  ranges  of 
velocity.  It  is  interesting  to  note  that  the  ratio  does  not  depend  upon  x0  /  p.  As  expected, 
because  of  thermal  coupling  of  mechanical  properties,  the  ratio  diminishes  when  velocity 
increases.  This  trend  is  directly  due  to  the  temperature  softening  of  the  flow  stress.  At  very 
high  velocities  the  ratio  reaches  very  low  values,  at  V  =  1000  m/s  the  asymptotic  value  is  less 
than  10%  in  comparison  to  the  ratio  at  V  =  1.0  m/s.  The  same  trend  is  shown  in  Fig.20  where 
the  ratio  of  the  coefficients  of  friction  is  shown  as  a  function  of  the  nominal  strain  rate 


t  =  V/h  in  the  logarithmic  scale.  The  asperity  size  h  =  20  pm  combined  with  velocities  V 
from  1.0  m/s  to  1000  m/s  yields  extremely  high  strain  rates,  from  104  1/s  to  107  1/s.  Those 
calculations  were  performed  for  different  fractional  area,  from  0.2  to  0.9  and  for  constant  non- 
dimensional  pressure  0.9. 
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A  general  definition  of  the  coefficient  of  friction  which  stems  from  the  asperity  model  is 
given  by  Eq.(19).  The  simplifications  discussed  previously  have  lead  to  the  definition  of  the 
coefficient  of  friction  based  on  plasticity  in  the  form  of  Eq.(23).  The  most  important  part  of 
the  analysis  is  the  question  how  to  introduce  the  ratio  x  /  p.  When  the  visco-plasticity  is 
introduced  in  Eq.(23)  in  the  form  of  the  constitutive  relations  derived  in  this  study,  Eqs  (34), 
(35)  and  (36),  the  ratio  x  /  p  becomes  strain-rate  and  temperature  dependent. 

Three  possibilities  based  on  plasticity  are  discussed  below. 

1.  According  to  the  Coulomb  hypothesis  the  tangential  force  (shear  stress)  increases 
proportionally  to  the  normal  force  (normal  stress).  Thus,  x  =  x  xy  and  p  =  xoy, 
0<x<l,  and  assuming  the  Huber-Mises  yield  condition  (isotropic  hardening), 

<7y  =^3ry,  the  ratio  x /pis  reduced  to  tip  =  1/^3  or  t!p  =  0.577 .  Thus,  the 

upper  limit  of  the  coefficient  of  Coulomb  friction  based  on  isotropic  plasticity  is  0.577. 

2.  Another  possibility  introduced  in  [1.4]  is  the  kinematic  hardening  where  the  yield 
surface  is  translated  in  the  direction  of  plastic  flow,  that  is  in  the  direction  of  shear. 

This  is  the  Prager-Ziegler  hardening.  Then,  the  ratio  x  /  p  =  z  /  <jy  ,  with  ay  =  const. 

3.  The  last  possibility  considered  in  this  study  is  an  assumption  that  the  pressure  p 
increases  proportionally  to  the  yield  stress,  p  =  x  ay ,  and  shear  stress  remains 

at  the  threshold  level  xo.  Again,  the  Huber-Mises  yield  condition  is  applied.  The  main 
physical  motivation  is  that  independently  of  the  level  of  the  normal  stress  the  micro 
-asperities  are  always  cut  by  adiabatic  shearing. 

Those  three  cases  lead  only  to  different  initial  levels  of  the  coefficient  of  friction  and 
different  effect  of  pressure,  but  they  do  not  change  the  general  trends  of  how  the  coefficient  of 
friction  evolves  when  the  velocity  of  sliding  increases. 

For  cases  1  and  3  the  equation  for  p,  (Eq.(49)),  is  reduced  to  the  following  forms 


and 


for  case  1 


ju  = 


for  case  3 
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where  (*)  is  the  multiplicative  term  after  To  /  p  in  Eq.(49).  When  the  following  variables  are 

introduced:  p  =  Xa/3  r0  and  p0  =  x0  -J3  t0  in  Eq.(49),  the  explicit  expression  for  p  (case  3)  is 
obtained 
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Both  expressions  for  p  ,  the  first  one  called  coefficient  of  friction  (Coulomb  sense),  and  the 
second,  called  friction  parameter,  Eq.(51),  have  been  evaluated  numerically.  The  solution  for 
©b  ,  that  is  Eq.(48),  was  introduced  into  Eq.(49)  or  Eq.(51)  in  order  to  find  p.  In  Fig  21  are 
shown  the  results  of  such  calculations  for  the  case  of  Coulomb  plasticity  where  the 
coefficient  of  friction  was  calculated  at  different  velocity  levels,  i.e.,:  1.0  m/s  <  V  <  1000  m/s, 
as  a  function  of  non-dimensional  pressure  x/x0  and  for  constant  fractional  area  0.9.  As  in  the 


quasi-static  case  this  coefficient  of  friction  increases  when  pressure  increases  as  well. 
However,  at  high  sliding  velocities  the  coefficient  of  friction  remains  almost  constant  at  a  very 
low  level,  independently  of  pressure.  The  most  important  results  are  shown  in  Fig.22  where 
the  coefficient  of  friction,  in  the  sense  of  Coulomb  plasticity,  is  shown  within  three  ranges  of 
sliding  velocity,  (a)  1.0  m/s  <  V  <  10  m/s,  (b)  10  m/s  <  V  <  100  m/s  and  (c)  100  m/s  <  V  < 
1000  m/s.  Calculations  were  performed  for  different  fractional  area  from  0.2  to  0.9  and  at 
constant  non-dimensional  pressure  0.9.  The  level  of  the  coefficient  of  friction  at  low  velocities 
is  relatively  low,  around  0.04  at  V  =  10  m/s.  A  substantial  decrease  in  the  coefficient  of 
friction  is  observed  at  high  sliding  velocities  .  This  is  in  complete  agreement  with  the  physical 
intuition  and  literature  data.  At  velocity  V  =  1000  m/s  an  asymptotic  value  of  p  is  reached 
around  0.005.  Overall  the  low  values  of  the  coefficient  of  friction  is  due  to  multiplication  of 
the  particular  effects  in  the  normalized  form,  that  is  all  fractions  are  less  than  one  and  they  are 
multiplied. 

Another  series  of  calculations  have  been  performed  to  analyze  the  friction  parameter  (non 
Coulomb  friction)  defined  by  Eq.(51).  In  this  case  the  friction  parameter  is  decreasing 
function  of  non-dimensional  pressure.  This  is  shown  in  Fig.23  where  the  friction  parameter  is 
drawn  at  different  velocities,  1.0  m/s  <  V  <  1000  m/s,  and  at  constant  fractional  area  0.9.  The 


result  that  the  friction  decreases  with  pressure  is  consistent  with  experimental  findings,  for 
example  [24].  The  following  series  of  calculations  was  similar  as  for  the  Coulomb  plasticity 
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series,  with  the  same  three  ranges  of  velocities  from  1.0  m/s  to  1000  m/s.  The  results  are 
shown  in  Fig.24,  (a)  1.0  m/s  <  V  <  10  m/s,  (b)  10  m/s  <  V  <  100  m/s  and  (c)  100  m/s  <  V  < 
1000  m/s.  Since  the  non-dimensional  pressure  was  assumed  in  both  cases  at  the  same  level  0.9 
and  assumed  values  of  the  fractional  area  were  the  same,  from  0.2  to  0.9,  changes  of  the 
friction  parameter  as  a  function  of  velocity  are  very  similar  for  both  cases.  Because  the 
absolute  values  of  the  coefficient  of  friction,  and  friction  parameter,  predicted  by  this  model 
are  relatively  low  and  their  evolution  as  a  function  of  velocity  seems  to  be  correct  the  only 
solution  at  this  stage  of  study  is  to  introduce  a  phenomenological  quasi-static  coefficient  of 
friction  po  determined  at  V  =  1.0  m/s,  or  assume  another  low  velocity  of  sliding.  Estimation 
of  the  adjusted  values  of  the  coefficient  of  friction  is  provided  for  all  definitions  of  x  /  p  by  the 
expression 

A  =  Ao(*)  (52) 

where  (*)  is  the  multiplicative  term  in  Eq.(49)  or  Eq.(5 1). 

In  conclusion,  the  trend  in  evolution  of  friction  as  a  function  of  velocity  depends  in  both 
cases  upon  constitutive  relations  coupled  with  thermodynamics  of  HAZ  defined  as  the  time 
evolution  of  ©b. 

11.  Discussion  and  conclusions 

After  analysis  presented  in  the  previous  chapters  of  this  study,  it  appears  that  one  of  the 
most  difficult  problems  in  the  mechanics  of  KE  projectiles  penetrating  geological  or 
cementitious  target  materials  is  a  proper  understanding  of  the  frictional  properties  between  the 
projectile  surface  and  the  target.  The  model  developed  in  this  study  is  limited  to  a  very  simple 
definition  of  the  state  of  the  projectile  surface  as  a  set  of  uniformly  distributed  micro- 
asperities  with  the  same  active  height  h  =  20  pm.  Of  course,,  more  general  statistical  models  of 
micro-asperities  of  different  heights  can  be  envisaged,  but  such  an  approach  is  at  present 
premature.  The  simplest  way  to  approach  such  a  difficult  problem  is  to  use  fundamental  laws 
of  physics.  In  this  case  a  physically  based  constitutive  relation  has  been  proposed.  The  main 
problem  in  this  study  was  lack  of  experimental  data  on  visco-plastic  behavior  of  steels  at  very 
high  strain  rates,  up  to  5*107  1/s  and  at  temperatures  higher  than  half  of  the  melting  point.  It 
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appears  that  this  area  of  testing  would  be  very  useful  in  providing  new  data  for  further 
applications. 

Although  the  heat  conduction  equation  has  not  been  solved  numerically  in  this  study,  a 
useful  approximation  of  the  evolution  of  temperature  in  the  bulk  material  was  used  after  [14]. 
This  approximation  permitted  finding  the  closed  form  solution  for  ©b.  Because  evolution  of 
the  bulk  temperature  could  be  estimated  it  also  permitted  finding  the  evolution  of  the 
coefficient  of  friction  as  a  function  of  the  sliding  velocity.  For  every  definition  of  the 
coefficient  of  friction  a  substantial  decrease  of  the  resistance  to  sliding  has  been  found.  This  is 
consistent  with  the  physical  interpretation  because  a  thin  melted  layer  of  metal  on  the 
projectile  surface  provides  lubrication.  An  open  question  remains  as  to  the  role  of  the 
hydrostatic  pressure.  The  experimental  data  reported  in  [24]  seem  to  indicate  that  the 
hydrostatic  pressure  reduces  friction.  In  Fig.25  are  shown  the  re-analyzed  results  of  the  pin- 
on-disk  experiments  reported  in  [24]  for  projectile  steel  on  gun  steel.  The  coefficient  of 
friction  was  determined  as  a  function  of  velocity  multiplied  by  pressure.  For  example,  the 
experimental  data  suggests  that  at  constant  velocity  the  coefficient  of  friction  will  diminish  as 
a  function  of  pressure  as  predicted  by  Eq.(51).  Thus,  an  advantage  of  the  model  presented  here 
is  that  the  hydrostatic  pressure  can  be  taken  into  consideration  in  a  correct  way.  The  physical 

limit  of  pressure  is  the  yield  stress,  oy  =  yfd  T0  ,  with  x0  =  520  MPa  the  maximum  pressure 

is  about  9  kbar.  Above  this  pressure  level  the  pressure  effect  in  the  model  becomes  empirical. 
Albeit  implementation  of  the  friction  model  in  the  form  of  Eq.(51)  into  numerical  codes  is 
rather  a  too  early  move,  preliminary  trials  should  be  recommended. 
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Figure  captions 


Fig.  1  Approximation  of  surface  topography  by  regular  array  of  columns  with  random 
array  of  heights. 

Fig.2  Top  view  of  the  surface  topography,  the  hot  spots  are  randomly  distributed. 

Fig.3  Elementary  cell  with  an  active  micro-asperity  which  is  cut  by  adiabatic  shear  band. 
Fig.4  Non-dimensional  number  of  contacts  as  a  function  of  non-dimensional  pressure,  Eq. 

.  (22). 

Fig.5  Quasi-static  coefficient  of  friction  (Coulomb  type)  as  a  function  of  non-dimensional 
pressure  for  different  fractional  area  of  elementary  cell,  Eq.(23). 

Fig.6  Quasi-static  friction  parameter  (non-Coulomb  type)  as  a  function  of  non-dimensional 
pressure,  Eq.(24). 

Fig.7  Logarithm  of  shear  strain  rate  versus  asperity  height  h,  velocity  of  sliding  1 .0  m/s. 
Fig.8  Logarithm  of  shear  strain  rate  versus  asperity  height  at  high  velocities,  200  m/s  <  V  < 
1600  m/s. 

Fig.9  .Theoretical  values  of  pseudo  viscosity  near  the  supersonic  speed  of  dislocations. 

Fig.  10  Changes  of  the  relative  pseudo  viscosity  as  a  function  of  homologous  temperature, 
Eq.(34). 

Fig.l  1  Failure  of  a  micro-asperity  by  ASB  in  elementary  cell. 

Fig.  12  Time  to  reach  specific  homologous  temperature  for  different  ranges  of  velocity, 

(a)  1.0  m/s  <  V  <  10.0  m/s,  (b)  10  m/s  <  V  <  100  m/s,  (c)  100  m/s  <  V  <  1000  m/s. 
Fig.  13  Time  to  reach  specific  homologous  temperature  as  a  function  of  velocity,  (a),  (b),  (c) 

the  same  velocity  ranges  as  in  Fig.  12. 

Fig.  14  Verification  of  the  failure  criterion  for  different  ranges  of  velocities,  (a),(b),(c)  -  the 
same  velocity  ranges  as  in  Fig.  12. 

Fig.  15  Heat  affected  zone  (HAZ)  for  two  different  time  scales. 

Fig.  16  Evolution  of  the  bulk  homologous  temperature  as  a  function  of  the  sliding  velocity  for 
temperature-independent  pseudo-viscosity. 

Fig.  17  Evolution  of  the  bulk  homologous  temperature  as  a  function  of  sliding  velocity, 
complete  temperature  coupling;  (a)  0  <  V  <  1 .0  m/s,  (b)  1 .0  m/s  <  V  <  1 0  m/s, 

(b)  10  m/s  <  V  <  100  m/s,  (d)  1 00  m/s  <  V  <  1 000  m/s. 
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Fig.  18  Evolution  of  the  bulk  homologous  temperature  as  a  function  of  the  fractional  area  of 
asperity  a  /  A,  (a)  1.0  m/s  <  V  <  10  m/s,  (b)  10  m/s  <  V  <  100  m/s,  (c)  100  m/s  <  V 

<  1000  m/s. 

Fig.  1 9  Ratio  of  the  coefficients  of  friction  for  three  ranges  of  velocities,  (a)  1 .0  m/s  <  V  < 

10  m/s,  (b)  10  m/s  <  V  <  100  m/s,  (c)  100  m/s  <  V  <  1000  m/s. 

Fig.20  Ratio  of  the  coefficients  of  friction  as  a  function  of  the  nominal  strain  rate  at  different 
fractional  area  and  constant  pressure. 

Fig.21  Coefficient  of  friction,  in  the  sense  of  Coulomb  plasticity,  as  a  function  of  the  non- 
dimensional  pressure  at  different  levels  of  velocity  and  constant  fractional  area  0.9, 
range  of  velocities:  1 .0  m/s  <  V  <  1000  m/s. 

Fig.22  Coefficient  of  friction,  in  the  sense  of  Coulomb  plasticity,  as  a  function  of  sliding 
velocity  V  for  different  levels  of  fractional  area  and  constant  non-dimensional 
pressure,  (a)  1.0  m/s  <  V  <  10  m/s,  (b)  10  m/s  <  V  <  100  m/s,  (c)  100  m/s  <  V  <  1000 
m/s. 

Fig.23  Friction  parameter  (non-Coulomb  friction)  as  a  function  of  the  non-dimensional 
pressure  for  different  levels  of  velocity  and  constant  fractional  area  0.9,  range  of 
velocities:  1.0  m/s  <  V  <  1000  m/s. 

Fig.24  Friction  parameter  (non-Coulomb  friction)  as  a  function  of  sliding  velocity  V  for 

different  levels  of  fractional  area  and  constant  non-dimensional  pressure,  (a)  1.0  m/s 

<  V  <  10  m/s,  (b)  10  m/s  <  V  <  100  m/s,  (c)  100  m/s  <  V  <  1000  m/s. 

Fig.25  Re-analyzed  data  for  friction  coefficient  reported  in  [24]  determined  for  projectile 

steel  on  gun  steel  by  high  speed  pin-on-disk  device;  range  of  pressure  0.24  kbar  <  p  < 
1.34  kbar,  range  of  velocity  3.0  m/s  <  V  <  550  m/s. 
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